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Abstract. We present several new sets of grids of model stellar atmospheres computed with modified versions 
of the AFLAS9 code. Each individual set consists of several grids of models with different metallicities ranging 
from [M/H] = -2.0 to +1.0 dex. Fhe grids range from 4000 to 10000 K in T cB and from 2.0 to 5.0 dex in 
loggr. Fhe individual sets differ from each other and from previous ones essentially in the physics used for the 
treatment of the convective energy transport, in the higher vertical resolution of the atmospheres and in a finer 
grid in the (T B s, logy) plane. Fhese improvements enable the computation of derivatives of color indices accurate 
enough for pulsation mode identification. In addition, we show that the chosen vertical resolution is necessary 
and sufficient for the purpose of stellar interior modelling. Fo explain the physical differences between the model 
grids we provide a description of the currently available modifications of AFLAS9 according to their treatment 
of convection. Our critical analysis of the dependence of the atmospheric structure and observable quantities on 
convection treatment, vertical resolution and metallicity reveals that spectroscopic and photometric observations 
are best represented when using an inefficient convection treatment. Fhis conclusion holds whatever convection 
formulation investigated here is used, i.e. MLF(a = 0.5), CM and CGM are equivalent. We also find that changing 
the convection treatment can lead to a change in the effective temperature estimated from Stromgren color indices 
from 200 to 400 K. 

Key words, stars: stellar atmospheres - stars: fundamental parameters - stars: 8 Scuti stars - physical data and 
processes: convection 



1. Introduction 

Convective transport of energy in a stellar atmosphere is 
one of the most complex astrophysical problems. Many 
of the approximations usually admitted for the stellar in- 
terior, such as diffusive radiative transfer, are no longer 
valid. Moreover, throughout most of a convective stellar 
atmosphere radiative losses are large enough to make con- 
vection less efficient in transporting energy than radiation. 
Only stars which have a surface convection zone (CZ) ex- 
tending deep into the stellar envelope can maintain effi- 
cient convective energy transfer near the bottom of their 
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atmosphere. On the other hand, inefficient convection ap- 
pears in all stars near the boundary of a convection zone 
close to locally stable regions. The modelling of inefficient 
convection requires a detailed knowledge about the effect 
of radiative gains and losses on the fluid flow. The sit- 
uation is particularly complex for stars which are cool 
enough to develop a granulation pattern, such as the sun. 
In this case, at identical geometrical depths, vastly differ- 
ent physical conditions may be encountered depending on 
whether upflow in a granule or downflow in an intergran- 
ular lane is considered. The former may be optically thick 
while the latter is already optically thin, a consequence 
of the extreme temperature sensitivity of the dominant 
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opacity source in the solar photosphere, the H ion (cf. 
also Stein & Nordhmd 1998). 

Currently, only very simple convection models are 
available for routine computation of extended grids of 
model atmospheres, while detailed numerical simulations 
are still unaffordable for applications that require the cal- 
culation of many thousands of individual model atmo- 
spheres over the HR diagram. 

Our intention here is first to review the convection 
models which are available for use together with the pop- 
ular ATLAS9 model atmosphere code by Kurucz (1993, 
1998) (see also Castelli et al. 1997). We provide an 
overview on what is known about the effects of the differ- 
ent convection treatments on model atmosphere structure 
and consequently on observable quantities. 

The second purpose of the present paper is to deter- 
mine to what extent the precision of fundamental param- 
eters derived from the observed stellar spectrum, i.e. T e g, 
gravity and metallicity depend on the model atmosphere. 

Another objective is to obtain very accurate colors and 
more importantly very accurate derivatives of colors, color 
indices and limb darkening coefficients. These quantities 
are needed in the procedure of pulsation mode identifi- 
cation which is the first and a crucial step in any seis- 
mological study. Indeed probing the stellar interior of a 
pulsating star requires the knowledge of the resonant cav- 
ity within which each mode propagates, i.e. the physical 
nature of the pulsation mode associated with each ob- 
served oscillation frequency. One such procedure is based 
on the computation of oscillation amplitude ratios and 
phase differences which in turn depend on the variation 
of the colors with effective temperature and gravity. The 
results of this application of the model atmosphere grids 
will be presented in the next papers of this scries (Barban 
et al. 2002; Garrido et al. 2002). Finally, due to their en- 
hanced resolution the new model grids are also useful to 
improve the outer boundary conditions of stellar struc- 
ture calculations (Montalban et al. 2001; DAntona et al. 
2002). 

These goals are part of a program performed in the 
framework of preparing the COROT space mission (see 
COROT web site). To achieve these purposes, we have 
used the ATLAS9 code in several versions modified for 
the convection zone treatment to compute new grids of 
model atmospheres, corresponding fluxes, surface intensi- 
ties, uvby colors, synthetic spectra for some representative 
lines, and compared them with relevant observations. We 
have three versions of the ATLAS9 code at our disposal: 

1. The original version from CDROM13 of Kurucz 
(Kurucz 1993) in which the convection zone is treated 
using mixing length theory (MLT). While ATLAS ver- 
sions from 5 to 8 remained basically close to the for- 
mulations given in Bohm-Vitense (1958) and in Cox 
& Giuli (1968), some improvements were added in 
ATLAS 9 (cf. Castelli et al. 1997). In Sect. 2 we dis- 
cuss the reasons for our specific selection among these 
improvements. 



2. The other two versions were provided by one of the 
authors (FK) who modified the code to include tur- 
bulent convection models from Canuto & Mazzitelli 
(1991, CM), and from Canuto, Goldman, & Mazzitelli 
(1996, CGM). 

Each convection model has been extensively used in the 
model atmosphere grid computations which we describe 
below. All the convection models are of local type and thus 
require the prescription of a characteristic length scale. 
Formally, it is possible to interchange the different length 
scales associated with the convection models. The motiva- 
tion for doing so and a particular example will be discussed 
in the next paper of this series (Kupka et al. 2002). 

This paper is organized as follows. In Sect. 2 we review 
previous works about the effect of the model structure on 
theoretical photometric colors and justify the need for new 
grids of model atmospheres. In Sect. 3 we describe the 
specific different convection treatments used and discuss 
their physical content. 

In Sect. 4 we give details of the grid computations. In 
Sect. 5 we set out and comment the role of the convec- 
tion treatments and convection parameters on the model 
structure, as well as its dependence on effective temper- 
ature, surface gravity, and metallicity. Finally, we discuss 
the consequences on observable quantities such as Balmer 
line profiles, flux distributions, and colors. 

2. A need for new grids 

The original grids of model atmospheres and colors based 
on the ATLAS9 code were published by Kurucz (1993). 
They were computed using the classical mixing length the- 
ory. Kurucz chose and fixed the mixing length parameter 
a, i.e. the ratio l/H p of convective scale length I and lo- 
cal pressure scale height H p , to be 1.25. He also used a 
prescription for overshooting at the top of the convection 
zone (cf. also Sect. 3) to achieve a better match between 
computed and observed solar fluxes for the range of a 
considered. The parameters obtained from the comparison 
with solar data were used for the entire grids published in 
Kurucz (1993). These grids have now been superseded by a 
new set with a slightly modified prescription of the over- 
shooting treatment (for details see Castelli et al. 1997). 
More recently, they have also become available in elec- 
tronic form (Kurucz 1998). 

Castelli et al. (1997) compared Johnson colors and the 
(b— y) and c indices from the Stromgren system with colors 
from grids of model atmospheres based on MLT with and 
without the overshooting prescription, and with an iden- 
tical choice for the mixing length. Considering different 
methods of determining T e g they concluded that models 
without the overshooting treatment yield more consistent 
results, while for the solar case a model with overshoot- 
ing was favored. As a consequence of this study, new grids 
of models, fluxes and colors were computed by Castelli 
without any overshooting for several metallicities and dif- 
ferent microturbulent velocities. They are available at the 
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Kurucz website ( "NOVER" grids). Castelli (1999) anal- 
ysed synthetic Johnson UBV colors from these model at- 
mosphere grids, all based on MLT with a = 1.25. She anal- 
ysed the effect of metallicity and microturbulent velocity 
and concluded that the indices are affected by both the 
convection treatment and the amount of line blanketing. 
This has to be considered in parameter determinations for 
stars with unknown metallicity. 

Kiinzli ct al. (1997) have used the revised version of 
model atmosphere grids of Kurucz (1998) to provide a 
new calibration of Geneva photometry for B to G type 
stars. Comparing their photometrically determined T e g 
and log g with evolutionary tracks for the Hyades they no- 
ticed a systematic trend in T e g below 7000 K and a rather 
pronounced "bump" in logg located in the same region. 
Both results were considered to indicate shortcomings in 
the model atmospheres used for the computations of the 
synthetic color indices. 

Smalley & Kupka (1997, SK) were the first to study 
the role of different convection treatments implemented in 
the ATLAS9 code of Kurucz (1993, 1998) for the synthetic 
uvby colors. They compared observed color indices with 
synthetic ones computed using two versions of ATLAS9: 
the original version of Kurucz (1993) based on MLT treat- 
ment of convection, with and without the overshooting op- 
tion, and another version modified to employ the convec- 
tion model of Canuto & Mazzitelli (1991, 1992), known as 
the CM model and described in Sect. 3. For the MLT they 
prove that models built with overshooting at the top of the 
convection zone, as illustrated in Castelli et al. (1997), 
are discrepant with the observed color indices. This con- 
firmed similar conclusions drawn by van 't Vcer-Mcnncret 
& Megessier (1996, hereafter VM) for the case of Balmcr 
line profiles. SK also showed that the CM models give 
results generally superior to those obtained with MLT us- 
ing a — 1.25, because they are in better overall agreement 
with the observed indices (6 — y)o and cq. The metallicity 
index nio was found to be the most discrepant one with ob- 
servations, the CM models remaining in good agreement 
only for stars with T e g larger than 7000 K, but clearly 
discrepant for solar type stars. A peculiar feature in the 
gravity sensitive co-index for T e g around 7000 K was found 
to be present in colors predicted using any of the convec- 
tion models investigated, similar to the results found by 
Kiinzli ct al. (1997) for MLT model atmospheres for the 
Geneva photometric system. 

A similar investigation to the one of SK for the 
Stromgren photometric system was done later by Schmidt 
(1999), but for the Geneva system. Moreover, he extended 
it to the CGM convection model which had meanwhile 
been implemented into the ATLAS9 code (see Sect. 3). His 
main conclusion, similar to the one of SK, can be summa- 
rized as follows: synthetic color indices are more sensitive 
to the scale length used than to the particular convection 
model. For instance, a value of a = 1.25 yields differences 
in the colors in comparison with models where a = 0.5 
which are much larger than the difference among CM and 
CGM models as well as MLT models with a = 0.5. He 



concluded that a value of a = 1.25 does not allow re- 
producing the observed photometric colors of late A and 
F stars. However, discrepancies were also found for the 
other convection treatments he had studied, in agreement 
with the results of SK on the uvby colors. 

Heiter et al. (1998) investigated the temperature struc- 
ture and observed quantities calculated with different con- 
vection models for two A Bootis stars with ([M/H], T e g) 
values of (-1, 6800 K) and (-2, 7800 K). They found a 
smaller difference between the synthetic colors and fluxes 
and the observations when using the CM model or MLT 
without overshooting compared to MLT with overshooting 
(a was set to 1.25 for the MLT models). For the cooler one 
among the two stars, the inclusion of overshooting changed 
the C, Ti, Cr, and Fe abundances derived from high res- 
olution spectra by +0.1 dex. They also compared the UV 
fluxes of these stars with IUE and TD1 measurements and 
found the CM convection model to yield results in best 
overall agreement while the discrepancies were largest for 
MLT models with a = 1.25 with overshooting. 

Recently, Gardiner et al. (1999) extended the compar- 
ison of SK to the CGM model for the case of Balmcr line 
profiles. It was found that differences between model at- 
mospheres based on the CM or CGM convection treat- 
ment, and models based on MLT without overshooting 
yield rather similar results, while MLT models with over- 
shooting are clearly different. A recommendation for a par- 
ticular model was found to be possible only for distinct, 
limited regions in T e g. Their results indicated that a more 
thorough study of the hydrogen line broadening mecha- 
nisms is necessary to draw more reliable conclusions on the 
convection model, as well as a larger number of standard 
stars with more accurately known fundamental parame- 
ters. For cool dwarf stars such as the sun, one source of 
problems in matching observed Balmer line profiles with 
synthetic ones has been to neglect the self-broadening (line 
broadening due to collisions with neutral hydrogen) in the 
hydrogen line profile calculations (Barklem et al. 2000). 
However, this effect is too weak in A and F stars to ex- 
plain the extent of the discrepancies found in matching 
the Balmer lines H Q and H^ with some of the model at- 
mospheres for the stars in the above mentioned works. 

From these previous works we have thus drawn the fol- 
lowing considerations for our grid computations. First, the 
overshooting prescription of ATLAS9 was generally found 
to be less successful in reproducing observations for A to 
G type stars, even though for solar observations the case 
is less settled. Thus, we have decided not to include mod- 
els computed with this treatment in our grids. However, 
for comparison we computed individual models with over- 
shooting (always using the correction by Castelli (1996)) 
for our case studies (Figs. 2, 5 and 8). 

Second, it has been found that model atmospheres 
which predict temperature gradients closer to the radia- 
tive one, i.e. where convection is less efficient than pre- 
dicted by MLT models with a > 1, arc in better overall 
agreement with observations. This was first noticed by 
Fuhrmann ct al. (1993) and, quite independently, for the 
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case of ATLAS9 models by VM where in order to repro- 
duce the sequence of Balmer line profiles of the sun with 
the same solar model they had to reduce the value of a 
of their MLT model atmospheres down to 0.5. Similar re- 
sults were found by Fuhrmann et al. (1993), VM and van 
't Veer-Menneret et al. (1998) for a large range of metal- 
licities and stars of spectral types between A5 and G5 
where Balmer lines are both strong and primarily sensi- 
tive to the temperature stratification. As shown above, 
this overall conclusion can also be drawn from other types 
of measurements such as photometry and is found to hold 
in particular for A type stars with T e g larger than 7000 K, 
while results for stars with lower T e g were generally more 
discrepant. Consequently, we have decided to base the ma- 
jority of our model grid computations on convection treat- 
ments which predict less efficient convection than the pre- 
vious model grids published by Kurucz (1993, 1998) and 
Castelli (1999). 

As far as oscillation mode identification procedures are 
concerned, it has been demonstrated that the dependency 
of the colors on T e ff and logg is not captured smoothly 
enough by the standard ATLAS9 models. The effects of 
the non smooth behavior of the color and limb darkening 
coefficient derivatives are larger than the expected effect 
used for identifying the modes (Garrido 2000). In order 
to obtain smooth variations of these quantities, we have 
found that it is necessary to compute our model atmo- 
spheres with a higher resolution in temperature distribu- 
tion with depth and built finer grids in T c g and \ogg. 

3. Convection Treatment: MLT and FST versions 
of ATLAS9 

3.1. Mixing length theory (MLT) 

Model atmospheres computed with ATLAS9 are based 
on the classical assumptions of stationarity and horizon- 
tal homogeneity. With these restrictions only some of the 
properties of stellar convection can be taken into account. 
ATLAS9 permits to include: 

1. the thermal contributions of convection to the energy 
flux through the atmosphere; 

2. the effect of convective motions on the line opacity due 
to the additional Doppler broadening of spectral lines 
caused by turbulent velocity fluctuations on length 
scales smaller than the mean free optical path. This is 
achieved by specifying a microturbulent velocity v m icro 
(cf. Gray 1992). 

3. Optionally, ATLAS9 permits to account for changes 
in pressure stratification due to a turbulent pressure 

Pturb- 

The convective energy flux F conv in the different versions 
of ATLAS has been computed traditionally with the clas- 
sical mixing length theory (cf. Biermann 1948; Opik 1950; 
Bohm-Vitense 1958; Castelli et al. 1997). Classical MLT 
includes radiative cooling of the fluid which is particu- 
larly important where convection is most inefficient, near 



the boundary of stably stratified layers. It requires the 
specification of a characteristic scale length I which is pre- 
scribed to be a fraction a of the local pressure scale height, 
P I 



H, 



(1) 



99 a 

I is used to describe the distance which fluid elements can 
travel before they dissolve. It also specifies the geometri- 
cal size of the fluid elements ( "bubbles" ) together with a 
second parameter, the ratio of the fluid element volume 
V over its surface area A. The quantity V/(Al) has been 
changed during upgrades of the ATLAS code (cf. Castelli 
1996). The present choice results in the same convective 
efficiency as the original one of Bohm-Vitense (1958) if 
slightly smaller values of a are used, i.e. the usual choice 
of a = 1.25 in the grids of Kurucz (1993, 1998) corre- 
sponds to an "aBv" of about 1.4 for A to G type main 
sequence stars. A detailed summary of the modifications of 
MLT as used in the ATLAS code can be found in Castelli 
(1996), together with various numerical coefficients which 
we have kept unaltered. 

One strong motivation to apply a more complete de- 
scription of stellar turbulent convection stems from the 
result that low values of the scale length parameter a, e.g. 
0.5, are required to fit Balmer line profiles for the sun and 
other cool dwarfs (Fuhrmann et al. 1993, VM), while much 
larger values (between 1 and 2) are necessary to reproduce 
their observed radii (Morel et al. 1994). Likewise, the scale 
length ratio has to be varied over an even larger domain 
(1 < a < 3) to reproduce the red giant branch in HR dia- 
grams of galactic open clusters and associations for stars 
with masses ranging from 1 Mq to 20 Mq (Stothers & 
Chin 1997, 1995). 

3.2. Full spectrum turbulence (FST) convection 
models 

An alternative to MLT which can address these problems 
was introduced by Canuto & Mazzitelli (1991, 1992) and 
is referred to as the CM convection model. An improved 
version was proposed by Canuto et al. (1996) which is 
known as the CGM formulation. A main intention behind 
both models was to improve the physical description of 
convection while keeping computational expenses as low 
as for MLT. Both models achieve this goal by providing a 
gradient (diffusion) approximation for the convective (en- 
thalpy) flux: 

Fco„v = K t p = K ia , d TH-\V - V ad )$(S), (2) 

where iVad = AacT 3 / (3np) is the radiative conductivity, 
$ = K t /K ra d is the ratio of turbulent to radiative con- 
ductivity, and 

= (^) ] =TF- 1 (V-V ad ) .3) 



dz 



dz ) a d 



is the superadiabatic gradient. The convective efficiency S 
is given by 

ga v (3l A v 



S = Ra • Pr 



X 



(4) 
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Ra and Pr are Rayleigh and Prandtl numbers of the con- 
vective flow, a v is the volume expansion coefficient, and 
the meaning of the other symbols is standard. We recall 
here that the thermometric conductivity \ is related to the 
radiative conductivity through K rad = c p p\ and that v is 
the kinematic viscosity. The quantity S is a useful measure 
of efficiency for flows which feature a very low Pr number, 
as occurs in stellar convection, and for which hence the 
detailed dependence on v can be neglected in parameter- 
isations. This is possible because viscous processes act on 
much longer timescales than radiation (t x = l 2 /x) and 
buoyancy (tb = (gctvP)^ 1 ^ 2 ) which in turn are respon- 
sible for the energy balance in stellar atmospheres and 
envelopes. Thus, the convective efficiency in a star can be 
characterized using only (i x /£b) 2 = S. The latter can eas- 
ily be related to an efficiency definition more common in 
astrophysics (Cox & Giuli 1968) that uses the quantity 



r = l((i + s)V3_i), 



(5) 



where 



E = 4A 2 (V-V ad ) 



81 



Q 1 l 2 c p p 2 nl 2 
12acT s 



2H, 



■,(6) 



and in which Q = TV^idV/d^p = 1 - (<9m^/dlnT) P 
is the variable average molecular weight. Using this nota- 
tion the MLT of B6hm-Vitense (1958) can be viewed as 
a phenomenologically derived prescription to compute <& 
which reads 



$MLT = 9 E _! ^ (1 + S) l/2_ i y^ 



(7) 



as mentioned by Canuto & Mazzitelli (1991) who pointed 
out that alternatively the MLT can be understood as a 
one-eddy approximation made for the spectrum E(k) of 
turbulent kinetic energy (see also Lesieur 1990). The latter 
describes how the kinetic energy of the velocity field gen- 
erated by convection is distributed among different spatial 
scales fc _1 . Canuto (1996) has shown how MLT underes- 
timates the convective flux in the high efficiency regime 
(S 1) while it overestimates F conv in the low efficiency 
regime (S <C 1). 

Both the CM and CGM convection models attempt 
to overcome the one-eddy approximation by using a tur- 
bulence model to compute the full spectrum E{k) of a 
turbulent convective flow for a given S, but keep the as- 
sumption of horizontal homogeneity and the Boussincsq 
approximation used in MLT. Hence, they are also referred 
to as full spectrum turbulence (FST) convection models. 

In the case of the CM convection model, the so-called 
eddy damped quasi-normal Markovian (EDQNM) model 
(Orszag 1977) of turbulence is used to compute &(S). This 
model provides a rather detailed treatment of the non- 
linear interactions in a turbulent flow, but requires the 
specification of a growth rate. The latter was computed 



from the linear unstable convective modes. To avoid the 
solution of the equations of the turbulence model each 
time in a stellar code, the results for F conv were tabulated 
in a dimensionless form. This was achieved by comput- 
ing the quantity $(Ra, Pr) for a large range of Ra and 
Pr numbers. For Pr < 10~ 3 the function <3? was found to 
saturate. This agrees with the previous remark that S is 
a useful measure of convective efficiency in a star, where 
Pr is even orders of magnitudes lower, and it was hence 
sufficient to consider only the results for the lowest Pr 
number for a tabulation of $(5), or actually $(£), given 
by the EDQNM model. Canuto & Mazzitelli (1991) found 
that can be represented by the following analytical 

fit formula to an accuracy of better than 3%: 



$ CM (S) = ai £ fc ((l + a 2 E) 



I)' 



where 



oi = 24.868, a 2 = 0.097666, 

k = 0.14972, m = 0.18931, n = 1.8503. 



(8) 



(9) 



The comparison with the CGM model published later is 
simplified if one considers a change of variable from S to 
S. In that case 

$ CM (S) =&iS fc ((l + & 2 S) m -l)", where (10) 
h = 14.288, b 2 = 0.0024115, 

and k, to, n are the same as above. (11) 

While the asymptotic behavior of both MLT and CM 
models are equal, i.e. they fulfill the limiting relations 
k + ran ~ 1/2 and thus 



$(S) ~ S 1/2 for S > 1 
as well as k + n ~ 2 and hence 
$(S) - S* 2 for S < 1, 



(12) 



(13) 



a distinguishing feature of <I> CM (S') is to yield about 10 
times more flux than (7) for S ^> 1, i.e. 



$ CM (S) - 10$ MLi (5) for S > 1 



while 



$ CM (S) ~ 0.1$ Mm '(5) 



MLT/ 



for 



S < 1. 



(14) 



(15) 



The function $ CM defined by (8)-(9) (or (10)-(11)) is 
only the first ingredient of the "CM model" . Because $ 
is computed as a function of local variables (4), it de- 
pends on a characteristic length scale which cannot be 
provided by the formalism itself. Following the physical 
argument that the Boussinesq approximation leaves no 
natural unit of length other than the distance to a bound- 
ary and that eddies near the boundary of the convection 
zone are smaller than in the middle of the same (stacking), 
Canuto & Mazzitelli (1991) proposed to take 



I = z 



(16) 
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where z is the distance to the nearest stable layer. The 
combination of (8)-(9) and (16) has subsequently been 
called the "CM model" . In this form it was implemented 
by Kupka (1996) into ATLAS9 and used for the model grid 
computations presented here, although other prescriptions 
of I had been implemented and experimented with as well. 

In a subsequent paper, Canuto et al. (1996) proposed a 
different FST convection model which avoided the usage of 
a growth rate. Rather, it was taken into account that the 
rate of energy input which feeds the velocity fluctuations 
and thus keeps convection from decaying is controlled by 
both the source of instability (buoyancy) and by the tur- 
bulence it generates. However, the treatment of the non- 
linear interactions had to be more simplified to keep the 
analytical model manageable. The equations of the turbu- 
lence model were solved in the limit for low Pr numbers. 
The new self-consistently computed input rate results in 
an increase of the convective flux for a given efficiency S 
which is largest at intermediate values of S ~ 300. For 
that reason a more complicated analytical fit formula had 
to be used to represent the predictions of the turbulence 
model to an accuracy better than 3% for all values of S. 
The CGM expression for $ reads 



$ CGM =F 1 {S)F 2 (S) 



(17) 



where Fi(S) has the same structural form as in the CM 
model, 

= (Ko / 1. 5f aS k ((l + bS) m -l) n , with (18) 



a = 10.8654, b = 0.00489073, 

k = 0.149888, m = 0.189238, n = 1.85011, 



while F2(S) is given by 
F 2 (S) = 1 



cSP eS r 
+ 



1 + dSi 1 + fS*' 



with 



(19) 



(20) 



c = 0.0108071, d = 0.00301208, 
e = 0.000334441, / = 0.000125, 
p = 0.72, q = 0.92, r = 1.2, t = 1.5. 



(21) 



Here, Ko is the Kolmogorov constant which has been taken 
1.7 in all our calculations, a value well inside of the ex- 
perimental range (Praskovsky & Oncley 1994). Note that 
(jjCGM ghowg the- same asymptotic behavior as $ MLT and 
$ CM in the limits of S < 1 and S > 1. Moreover, the 
CM and CGM functions <f> approach these limits in a very 
similar manner, because F2(S) — > 1 for both very large 
and very small S, and the power exponents k, m, n of (9), 
(11), and (19) are almost identical. However, while 



$ CGM (S)~$^(S) for S > 1, 



S CM 



(22) 



the low efficiency results differ, as 

$ CGM (5) ~0.3$ MLT (S) for 5 « 1 (23) 



(cf. (15)). On a logarithmic scale, the low efficiency limit 
of (17) is almost exactly the average of the fluxes of (7) 
and (10)— (11). The second difference between the two FST 
convection models is the choice of the scale length I which 
Canuto et al. (1996) have proposed to be 

l = z + a*H pMp . (24) 

This accounts for the observed fact that convection pen- 
etrates into neighboring stable regions and thus the 
scale length cannot decay to zero right at the layer 
where the stratification becomes stable according to the 
Schwarzschild criterion. The additional term in (24) is 
thus supposed to account for overshooting and provides 
a possibility for small adjustments, if exact stellar radii 
are needed, e.g. in helioseismology. However, the meaning 
of overshooting in this context must not be confused with 
the overshooting option offered by the ATLAS9 code. This 
point deserves special attention to which we turn in the 
following. 

3.3. Length scale parameters and overshooting 

The term a*H Pitop in (24) accounts for the increase of the 
efficiency of convection due to convective penetration at 
the boundary between a stably and an unstably strati- 
fied region compared to a rigid boundary, for instance a 
fixed plate. The stellar scenario thus implies to increase 
the scale length I which can no longer be forced to zero as 
in (16). The total flux within convectively stable layers is 
still taken equal to the radiative flux. On the other hand, 
the overshooting prescription included in Kurucz (1993, 
1998) as illustrated in Castelli et al. (1997) was invented 
to take into account that overshooting directly changes 
the temperature gradient also in a stable region next to 
a convection zone. The procedure suggested is to simply 
smooth out the convective flux over as much as 0.5 H p in 
each direction around the last point where V = V ra d- This 
mimics the well-known property found in many numerical 
simulations (e.g. Hurlburt et al. 1986, 1994) and in so- 
lutions of the nonlocal Reynolds stress equations (Kupka 
1999; Kupka & Montgomery 2002) where F conv > even 
though V — V a d < in layers right next to a neighboring 
convection zone. A steeply decaying F conv cannot be mod- 
eled this way while the adjacent region where -Fconv < 
has to be neglected by taking V = V ra d- The effect of 
this flux smoothing procedure of ATLAS9 on the emer- 
gent flux is large enough to provide an additional degree 
of freedom to improve the match of solar observations by 
adjusting the smoothing width. 

In the CGM model, the parameter a* of (24) is typi- 
cally of order 0.1 and may be slightly changed to compen- 
sate for uncertainties in opacities and in the treatment 
of convection. Values of 0.08 and 0.09, similar to Canuto 
et al. (1996), were used for the different grids presented 
in Sect. 4. However, the effect of such small changes is 
minute. No inconsistencies were found in a recent work by 
Montalban et al. (2001) when model atmospheres com- 
puted with a* =0.09 were matched on top of stellar en- 
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velopcs at different tr oss , despite a slightly larger value 
was used in the stellar structure computations to obtain 
the correct solar radius when using the most recent opac- 
ity data. On the other hand, using a* to compensate for 
the Boussinesq approximation and various homogeneity 
assumptions in ATLAS9 by a match to, say, the entropy 
jump near the stellar surface as found from numerical sim- 
ulations (cf. Ludwig et al. 1999 who used a combination 
of the CM fluxes (8)-(9) and the scale length (24)) may 
require larger variations for models very different from the 
sun. However, such a procedure cannot bring the temper- 
ature gradient of ATLAS9 model atmospheres into agree- 
ment with the simulations. The latter avoid horizontal ho- 
mogeneity assumptions but cannot be afforded together 
with a treatment of frequency dependent radiative trans- 
fer which is comparably sophisticated as that one used in 
ATLAS9. Hence, emergent fluxes, spectra, and photomet- 
ric colors will be different as well. As long as such a match- 
ing procedure is not shown to allow an improved match of 
fundamental star data over extended parts of the HR dia- 
gram (and thus improving over present models, cf. various 
publications discussed in Sect. 2), its practical advantages 
appear more limited. For that reason, we have preferred to 
use the CGM model as intended by its authors and stud- 
ied grids with a constant a* which makes them suitable 
to be matched with stellar structure calculations using the 
same treatment of convection (Montalban et al. 2001). 

3.4. Implementation of FST models into ATLAS9 

In the ATLAS9 implementation of the CGM convection 
model the quantities (17)-(21) are actually computed as 
functions of S. Thus, only minimal changes were neces- 
sary in the subroutine TCORR, which performs the tem- 
perature correction, and in CONVEC, which computes 
the convective flux, for replacing the CM with the CGM 
model. TCORR and CONVEC were also the only subrou- 
tines that had to be changed for implementing the CM 
model into ATLAS9. The scale length of the CGM model 
is evaluated in the following way: 

I = min(z t0 p + a*H p , tov , z bo ttom + a* Hp, .bottom) (25) 

This choice makes convection slightly more efficient in 
comparison with (24) and more consistent with the idea of 
accounting for overshooting, as the latter is also expected 
to occur below convection zones. For most model atmo- 
spheres we found that the differences between these alter- 
native prescriptions are either zero or negligibly small, be- 
cause the temperature gradient for convection zones which 
are entirely contained within the atmosphere is practically 
radiative while for convection zones extending below the 
atmosphere the evaluation of I in. a pure model atmosphere 
code necessarily has to occur at the top of the convection 
zone. 

We note here that in principle (8)-(9) and (17)-(21) 
could also be used together with the common scale length 
I = aH p with a < 1, or other scale lengths. Results on 
such calculations will be reported in Kupka et al. (2002). 



3.5. Turbulent pressure and the optically thin limit 

For the CM model, a prescription for the turbulent pres- 
sure was published as well, although the results were given 
only for S ^> 1 and in tabular form. In stellar atmo- 
spheres, S 3> 1 is usually attained only in cool stars and 
close to the bottom where the Rosseland mean optical 
depth tross > 10. Hence, the ATLAS9 implementation of 
the CM model does not account for turbulent pressure. 
On the other hand, for the CGM convection model ana- 
lytical fit formulae for v tur b arid pturb were published by 
Canuto et al. (1996) which can be used even for S -C 1 
and were implemented into ATLAS9 as well. A number 
of model atmospheres for A to early M type dwarfs and 
for giants were computed with the CGM model with and 
without the prescription of p t urb- Differences were found 
only for stars with deep envelope convection zones, al- 
though in most cases both T and P changed by less than 
0.1% for tr oss < 5, and by no more than 0.5% to 1% for 
10 < trqss < 100. As the inclusion of pturb slowed down 
the convergence of models while spectra and colors re- 
mained indistinguishable from the case pturb = 0, all the 
CM and CGM model atmospheres grids presented here 
are computed without a pturb, just as their MLT counter- 
parts. We note that for stellar structure calculations the 
change in temperature structure due to pturb may be more 
important than for flux predictions derived from ATLAS9 
model atmospheres. To avoid discrepancies with the CGM 
treatment as used in the model grids a reasonable compro- 
mise is to match model atmospheres and stellar envelopes 
at a r Ross — 10. 

Following a suggestion by Canuto (private communi- 
cation) the correction of Spiegel (1957) for radiative losses 
in optically thin media was implemented for the case of 
the CM model. However, except for late K and early M 
dwarfs, where ATLAS9 models are not reliable any more 
due to the dominance of molecular lines, the effects were 
found to be negligible. The primary reason for this are 
the very low values of -F CO nv predicted by the CM model 
for t Ross < 2 for stars with T off > 4000 K. For the CGM 
model, convection is slightly more efficient, but still the 
effects of such a correction arc expected to be very small. 
Therefore, no further experiments with radiative loss rates 
were made with FST convection models. The case is dif- 
ferent for MLT where the results are more sensitive to the 
different cooling rates of "optically thin bubbles" , as 

$ MLT (S) > $ CGM {S) > <S> CM (S) for S < 10 (26) 

because of (15) and (23) and due to the much larger / 
of MLT for a > 1 if z < H p . A correction of F conv for 
the optically thin limit is always included in the MLT 
implementation of ATLAS9 (Kurucz 1993; Castelli et al. 
1997). 

4. Model grid computation 

Two model grids have been computed independently at 
the Paris and Vienna observatories. 
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At the Paris Observatory an automatic procedure was 
created by one of the authors (DK). The procedure is in- 
teractive, and allows the computation of grids of model 
atmospheres based on the ATLAS9 code, of Balmer line 
profiles, surface fluxes and intensities, colors and synthetic 
spectra, all in one run. The flux and temperature compu- 
tations are iterated until the following convergence criteria 
are satisfied: the maximum of the flux and flux derivative 
errors have to be equal to or less than one and ten percent, 
respectively. In addition, the maximum of the temperature 
correction has to be equal to or less than one K. 

In the MLT case, we started from the original Kurucz 
grids (Kurucz 1993, 1998) and recomputed the models by 
the scaling procedure of the ATLAS9 code. The thickness 
of the layers of the model atmospheres was divided by 2 
or 4 in comparison with the original Kurucz (1993, 1998) 
models, in order to solve numerical instabilities in the iter- 
ation procedure for the flux computation, and to provide 
more accurate photometric colors (see Sect. 5.2.2 and next 
paper in this series). Models with higher resolution con- 
verged faster and smaller flux errors were achieved. 

The parameters used for these model grids are given 
in Table 1. We recall that the metallicity is given in 
terms of the logarithmic ratio between the total number 
of atoms of each species, except for hydrogen and helium, 
over the number of hydrogen atoms, with respect to the 
solar metallicity defined in the same way. For instance, 
[M/H]=0.0 and —1.0 means that the opacities entering 
the model calculations are computed using either solar el- 
ement abundances or solar element abundances divided 
by 10 for all elements other than hydrogen and helium. 
The MLT models were computed for two values of a, the 
original value used by Kurucz a = 1.25, and the lower 
value a = 0.5, chosen for reasons given in Sects. 3 and 5. 

In the CM and CGM cases, we started from our MLT 
models with a = 0.5, and computed grids with the same 
set of parameters. For the CGM convection a value of 
a* = 0.08 was chosen (see Sect. 3 for a discussion). 

At the Vienna Observatory, model grids with several 
combinations of convection treatment and vertical resolu- 
tion were computed for slightly smaller step sizes in T e s, 
larger step sizes in logg and more [M/H] values. For MLT 
models a value of 0.5 has been chosen for a. Convection 
has been turned off for models with T e g >8600 K, because 
the convective flux can be neglected for higher temper- 
atures, as can be seen from Fig. 4. As in the Paris grid 
the uppermost layer is located at log tr. oss = —6.875. The 
difference of consecutive layers in log TRossis 0.125 and 
0.03125 for models with 72 and 288 layers, respectively. 
In addition to the model atmospheres, fluxes and colors 
in 12 systems have been computed. Furthermore, infor- 
mation on the convergence extracted from the ATLAS9 
output is provided for each model. The atmospheric and 
computational parameters arc summarized in Table 1. 

The grid computations were performed with the perl 
package SMGT (Stellar Model Grid Tool), described 



in Schmidt (1999) 1 . This non-interactive program runs 
ATLAS9 repeatedly until the convergence criteria are sat- 
isfied for each model. The output of ATLAS9 is evalu- 
ated directly and selected information is provided for each 
model, such as the root mean square (RMS) and maximum 
values of the flux and flux derivative errors, the maximum 
of the convective to total flux ratio, the extension of the 
convection zone, and the optical depth where the temper- 
ature equals T e gf. The grids defined in Table 1 are available 
on CDROM on request from the authors. 

We note here that two different, but overlapping grids 
of model atmospheres were computed as there were dif- 
ferent applications in mind. The main motivation for the 
computation of the Paris grids was to calculate photomet- 
ric colors and their derivatives with respect to T e gf and 
log g, which will be used in view of seismic applications 
(Watson 1988; Garrido ct al. 1990; Balona & Evers 1999). 
This required rather small steps in logg, but a restricted 
range for T e g and few metallicity values. The results of this 
specific application will be discussed in a subsequent pa- 
per of this series. The Vienna grids, on the other hand, are 
intended for general use, which is the reason for choosing 
intermediate values for the parameter step sizes and cover- 
ing as much of the HR diagram as possible. Examples for 
already published applications of these grids can be found 
in Montalban et al. (2001, see below) and in D'Antona 
ct al. (2002). 

4.1. Resolution 

To show that for specific applications it is necessary to 
use the models with 288 layers, we examined the quan- 
tity Az = z(tr oss = 3.162) — ^(i^conv = 0), where z is the 
depth (distance from top layer) in the atmosphere in km 
and z(F conv = 0) is the depth of the upper limit of the con- 
vection zone. This quantity has been used by Montalban 
et al. (2001) for the calculation of the convective scale 
length in stellar interior models which use convective at- 
mospheres computed with ATLAS9 as a boundary condi- 
tion. It turned out that for a particular region in the HR 
diagram, calculating this quantity from atmospheric mod- 
els with 72 layers results in unphysical oscillations in solar 
evolutionary tracks which disappear for higher resolutions 
(J. Montalban, private communication). Fig. 1 shows the 
values of Az for a small grid of CGM model atmospheres 
with [M/H]=0, T off =4200. . . 4800 and log #=3.0. . . 3.6 for 
four different resolutions, with the stepsize A log tr oss di- 
vided by two for each successive resolution value. For 144 
layers, the results are rather different from the 72 layer 
ones (note the peak at (T off , logg) = (4400, 3.2)). There 
is a small change when increasing to 288 layers, whereas 
the change when using 576 layers is negligible. This shows 
that 288 layers are sufficient in low to moderate temper- 
ature atmospheric models, in particular as all structural 

1 A summary is given in the Appendix and direc- 
tions for the use of this program can be found at 
http://ams.astro.univie.ac.at/"" heiter/smgt_usage_l.html. 
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Table 1. Atmospheric and computational parameters of the model atmosphere grids. 





Paris 




Vienna 






Min Max Step 




Min Max 


Step 


T cft [K] 


6000 8500 250 




4000 10000 


200 


l°g 9 


2.0 4.5 0.1 




2.0 5.0 


0.2 


[M/H] 


-1.0, 0.0, +1.0 


-2.0, 


-1.5, -1.0, -0.5, 


-0.3, -0.2, -0.1, 






0.0, +0.1, +0.2, +0.3, +0.5, +1.0 


V micro [kmS _1 ] 


2 




b , l b , 2, 


4 


Convection 


MLT CGM CM 


MLT 


CGM CGM 


CM 


Parameter 


1.25, 0.5 0.08 


0.5 


0.09 




Alog TRoss a 


0.0625 or 0.03125 


0.125 


0.125 0.03125 


0.03125 


Number of layers 


143 or 285 


72 


72 288 


288 


a log TRoss(tOp) = 


-6.875 


b in preparation 






Fig. 1. Az = z(t Ro bb= 3.162) -z(F conv = 0) for a CGM model 
grid with four different numbers of layers equally distributed 
between log tr, ss =— 6.875 and +2. This quantity is used for 
the calculation of the convective scale length in stellar interior 
models and the graph shows its sensitivity to depth resolution 
in this part of the HR diagram. 

quantities (e.g. the temperature gradient V) are resolved. 
For models with T e g> 10000 K, on the other hand, we 
verified that 72 layers are sufficient. 

5. The effects of Convection Treatment 

5.1. Effects on model atmosphere structure 

We first examine changes of temperature and convective 
flux distribution when using different convection models. 
Figs. 2, 3, and 4 show the intricate dependence of the effect 
of convection treatment on T c g, logg, and [M/H] of the 
model. 

Fig. 2 displays the temperature and the convective flux 
as a function of Rosseland optical depth (log tr ss) cor- 
responding to the models used for three specific main se- 
quence solar metallicity stars which have been chosen so 
as to cover the temperature range of interest: the Sun, 
Procyon - a well studied reference star, and [3 Ari - a 
well observed hot star. For each star, several models are 
computed which differ only for the convection treatment. 
T e ff, logg, metallicity, and microturbulcnt velocity of the 
models are taken from previous detailed analyses (by CV 



for (3 Ari, van 't Veer-Menneret et al. (1998) for Procyon 
and Kurucz (1998) for the Sun). 

The slope of the T — t relation within the CZ indicates 
the efficiency of the convection transport. It is steeper for a 
less efficient convection, i.e. a temperature gradient closer 
to the radiative one. For instance, in Fig. 2a, it can be 
seen that the CM model predicts the least efficient con- 
vection, followed with increasing convective efficiency by 
the MLT (a = 0.5), the CGM and the MLT (a = 1.25) 
models. The same trend is observed for the convective flux 
in Fig. 2b. This is a consequence of the fact that radiative 
losses of the convective fluid arc always large within the 
stellar atmosphere where the gas is optically transparent. 
Hence, the inequality chain (26) always holds at least at 
lower optical depths (r < 1). The scale lengths (1), (24), 
(16) of MLT, CGM and CM respectively also obey such an 
inequality chain for distances z closer to stably stratified 
layers than aH p and for the ranges of a and a* consid- 
ered in our work. Therefore, the amount of convective flux 
and the associated T — t relations shown in Fig. 2 are an 
immediate consequence of these inequality chains. 

For hotter stars (T e ff~ 8000 K), as the convective effi- 
ciency decreases, all convection models predict a temper- 
ature gradient close to the radiative one. 

The effect of convection treatment on the atmosphere 
structure depends on metallicity, gravity and T e g in a com- 
plex way, as illustrated in Fig. 3. Evidently, a metal rich 
atmosphere reduces the efficiency of the convection trans- 
port, as does a low gravity, or a high T e fj. The influence of 
log g on the convective efficiency depends strongly on T e g, 
[M/H], and the convection model. For instance, a model 
at T c g = 6500 K, logg = 2.5, cumlt = 0.5, and ten times 
solar metallicity is completely radiative (see Fig. 3a, thin 
dashed lower line), while for logy = 4.5 and identical pa- 
rameters otherwise a small deviation from radiative strat- 
ification is found (thick dashed lower line). This deviation 
grows significantly when decreasing the metallicity to one 
tenth of the solar one (Fig. 3b). 

The variation of the maxima of the convective flux with 
T cff , log g, and [M/H] is shown in Fig. 4 for the CGM 
models. The decrease of convective flux with increasing 
metallicity is a consequence of the lower mass density (p) 
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log t Ross log t Ross 

Fig. 2. Distributions of temperature (left panels) and ratio of convective to total flux (right panels) with Rosseland optical 
depth for the three model atmospheres adopted for the Sun (a, b), Procyon (c, d) and the A5V star (3 Ari (e, f). For each star 
we show five models computed using different convection treatments. 




Fig. 3. Temperature versus Rosseland optical depth for models with two different log g values, computed with MLT without 
overshooting (a: [M/H] = +1, b: [M/H] = — 1, two different values of a for each [M/H]) and with FST convection formulation 
(c: [M/H] = +1, d: [M/H] = -1, CM and CGM for each [M/H]). The models are represented for two different values of T cff in 
each panel, using the same line styles: 7500 K for the upper four curves and 6500 K for the lower four curves. 



found in metal rich atmospheres. The latter is a result 
of the increased opacity, which requires a smaller column 
density for a given optical depth. Due to increased line 
blanketing in metal rich atmospheres, the requirements 
of flux constancy and hydrostatic equilibrium then result 
in both lower temperature and lower pressure in the out- 
ermost layers. As a consequence, lower densities are also 
found near the boundary of the CZ. This makes convec- 
tion less efficient, although this effect is partially compen- 



sated by a higher convective velocity found for metal rich 
atmospheres. Fig. 4 also shows the influence of using a 
four times higher resolution in optical depth, resulting in 
a much smoother run of the curves and in a small shift 
towards lower temperatures. 
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4000 4500 5000 5500 6000 6500 7000 7500 8000 8500 



Teff [K] 

Fig. 4. Maximum of convective to total flux ratios as a function 
of T e ff, for the following values of logg: 2, 3, 4, 5, represented 
by increasing line thickness, and three different metallicities as 
labeled for logg=2, with the same trends for the other values. 
The dashed curve was computed with [M/H]=0 and loggi=4, 
but with a four times higher depth resolution. The CGM con- 
vection model was used in all cases. 

5.2. Consequences on observable quantities 
5.2.1. Balmer line profiles 

The effect of changing the physical parameters entering 
the models on the Balmer line profiles (BLPs) is very 
complex. This is illustrated in Fig. 5, where the synthetic 
profiles for several different convection models are com- 
pared to the observed ones for the same three stars as 
in Fig. 2. The spectra shown in Fig. 5 were obtained at 
the Haute-Provence Observatory, with the spectrograph 
Aurelie attached to the 152cm reflector, equipped with a 
CCD receptor. The resolution is about 25000. The Aurelie 
spectra are observed in the first or second order, depend- 
ing on the wavelength. The wavelength range is 200 A, 
and the continuum tracing is local, using the most suit- 
able windows. With a signal to noise ratio of at least 400 
we can expect an accuracy for the continuum location of 
0.3 %, i.e. 0.5 % for the ratio of line to continuum fluxes 
in the line wings. This corresponds to a 30 to 60 K change 
in effective temperature for F to G stars. 

In the case of H a the effects are never larger than 0.5%. 
Fig. 5 shows that the H a profile is insensitive to the choice 
of any of the scale lengths or convection models discussed 
in Sect. 3, while in the case of Hp the profiles computed 
with MLT and a = 1.25 are too narrow. As an example, in 
the case of Procyon this Up profile must be computed with 
T e ff around 300 K higher to represent the observed profile. 
The insensitivity of H Q to any convection treatment is one 
of the reasons why it is a very good temperature indicator. 
However, it is formed close to the boundary of convectively 
instable layers and therefore can be modified by inclusion 
of overshooting. 

Fig. 5 is a convincing illustration that by the use of the 
CM or CGM convection treatment, the observed H Q and 



Hp profiles can be represented by the same atmosphere 
model, and this constraint can be achieved in the MLT 
case provided a value for a of about 0.5 is adopted. Thus, 
we want to emphasize that for the three stars investigated 
here, less efficient convection within ATLAS9 type model 
atmospheres allows the best fit of H a and Hp using the 
same atmosphere model. 

In the MLT case, the consequences of these effects on 
the BLPs have been extensively described by VM and 
Fuhrmann ct al. (1993), who demonstrated that the BLPs 
are sensitive probes of the atmosphere structure and of 
effective temperature for late A, F, and G dwarf stars. 
Fig. 6a, b illustrates the effect of changing the MLT pa- 
rameter a on the BLPs, and to what extent this modi- 
fication depends on the model parameters. This sensitiv- 
ity strongly depends on the selected combinations of the 
model parameters. For instance, the largest differences are 
seen for models with T e g between 7000K and 8000K, \ogg 
= 4.5 and high mctallicity. In contrast, the differences arc 
insignificant at low gravity, high metallicity and high tem- 
perature. The shape of the profiles is also affected, the 
most for low temperature and low metallicity models. 

Fig. 6c, d shows the difference between two line pro- 
files, computed with MLT and CM. The difference be- 
tween the CGM and CM LL3 profiles was not plotted, be- 
cause it is similar to MLT(a = 0.5) — CM. The differences 
with MLT(a = 1.25) are the largest and strongly depend 
on temperature and gravity, but less on metallicity. These 
statements mean that the most efficient convection treat- 
ment is MLT with a = 1.25, in agreement with the T — t 
laws shown in Figs. 2 and 3. From the observer's point of 
view, Fig. 6a-d also reveals that Balmer line profiles have 
to be measured and normalized to an accuracy of at least 
0.5% to draw a clear distinction between convection mod- 
els with different efficiency. Insufficiently determined pro- 
files may thus easily introduce erroneous trends or large 
scatter when analyzing their dependence on a particular 
convection treatment. 

We stress here that the sensitivity to MLT's parame- 
ter a strongly depends on gravity. For instance, for T e g 
< 7500 K all convection models with low gravity yield in- 
efficient convection. This is due to the fact that low grav- 
ity implies lower densities, and the convective efficiency 
is related to the density as explained in subsection 5.1 
above. Moreover, we suggest to consider the commonly as- 
sumed insensitivity of BLPs to gravity change for T e g be- 
low 8000 K with real caution. Indeed, Fig. 7 shows clearly 
that the sensitivity of BLPs to gravity changes depends 
more than usually expected on mctallicity, T cB , and fi- 
nally on all parameters playing a role in the efficiency of 
the convective transport. It depends also on the gravity 
itself, the second derivative is not zero. This effect is most 
important for the highest metallicity and largest T e ff- The 
main reason is that an increase of metallicity leads to a 
lowering of the density on the one hand while a higher ef- 
fective temperature favors radiative transfer on the other. 




Fig. 5. H a (left panels) and H/3 (right panels) line profiles computed with different convection models represented by the 
line styles as in Fig. 2, together with the observed profiles for the same three stars as in Fig. 2. 
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Fig. 6. Differences between normalized fluxes in Balmer line profiles obtained from two models differing only by the convection 
treatment, versus the distance from the line center AA in Angstrom units. Fluxes are normalized to the local continuum flux 
F c . The two left panels display the difference of fluxes computed using MLT models differing by the a parameter: F(a — 
1.25) — F(a — 0.5) for two metallicities. The different line styles correspond to different T fr's: full lines for 6500 K, dotted lines 
7500 K and short-dashed ones 8500 K, and thin and thick lines respectively correspond to log g — 2.5 and 4.5. In the right panels 
the models differ by the convection formulation, and are identifiable as follows: for MLT(q = 1.25)— CM full lines correspond 
to T eff = 6500 K, long dashed lines to T cB = 7500 K; for MLT(a = 0.50)- CM dotted lines correspond to T eS = 6500 K, and 
short-dash-dotted lines to T ef r = 7500 K. Thin and thick lines have the same meaning. 



5.2.2. Fluxes and colors 

We have computed fluxes for solar models with differ- 
ent convection treatments as follows: CGM, CM, MLT 
(a=0.5), MLT (a=1.25), MLT (a=1.25 with overshoot- 
ing) . The same parameters have been chosen for all mod- 
els: T off = 5777 K, log 5 = 4.4377, w micro = 1.5 kms" 1 , 
element abundances from Anders & Grevesse (1989), ex- 
cept for Fe, for which the current value ol\og(N-p c / N tot ) = 
—4.51 was used (Kurucz 1998). To compare the calcu- 
lated solar fluxes to observations, solar irradiance data 
have been taken from Neckel & Labs (1984), and in addi- 



tion from two more recent sources: Lockwood et al. (1992, 
Lowell Observatory, 1985) and Thuillier et al. (1998, 
SOLSPEC spectrometer on ATLAS I mission, 1992). The 
irradiances in the region of maximum emitted radiation, 
i.e. 410 to 510 nm, are displayed in Fig. 8. 

As can be seen, the three observational data sets are 
different from each other by up to 15 % (upper panel), al- 
though Neckel & Labs (1984) estimate 0.5 % as an upper 
limit for the local systematic error of their measurements. 
But they point out that intrinsic intensity variations de- 
pending on solar activity can occur when comparing mea- 
surements made at different times. These amount to 2 % 
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Fig. 7. Differences between Hp line profiles computed with 
log g — 2.5 and 4.5. Different convection models are represented 
by different line styles, for metallicities of —1 and +1 by thin 
and thick lines respectively, and for three effective tempera- 
tures. 



in certain spectral regions (e.g. the Call K line) in their 
data, which are derived from observations made over a 
20 yr period (see also Livingston ct al. 1991). For compar- 
ison, Lockwood et al. (1992) give an upper limit for the 
errors of their measurements of 2 % (their observations 



were made at a phase of low solar activity), and Thuillicr 
et al. (1998) quote a mean uncertainty of 2-3 % (data ob- 
tained at high solar activity) 2 . Detailed discussions of the 
error sources and comparisons with previous observations 
are given in each of the three references. 

However, the mean of the maximum relative difference 
between the irradiances from the three sources is 5 % in 
the region of 450 to 480 nm, which is much larger than 
the differences between the fluxes calculated with different 
convection models (2 %, cf. lower panel of Fig. 8). The CM 
and MLT (a=0.5) fluxes are almost identical to the CGM 
flux. Therefore, the solar irradiance measurements cannot 
be used to decide between the various models. A similar 
conclusion would result if measurements of solar central 
intensity would be used, because the error estimates by 
Neckel & Labs (1984) for these measurements (the other 
two sources did not include this kind of measurements) are 
equal to that for the irradiance spectra. Thus, we regard 
tests of central intensity calculations against observations 
(e.g. Castelli et al. 1997) as having limited significance 
until accuracy and absolute calibration of these data will 
have been established with the necessary reliability. 

The general dependence of the calculated flux on the 
convection model can be seen in Fig. 9, where the ratios 
of MLT and CGM to CM fluxes are displayed for two 
different values of T c s and logg and the extreme case of 
[M/H]= -1. For all cases, the CGM flux is closest to the 
CM flux, followed by MLT (a=0.5) and with a larger dis- 
crepancy by MLT (a=1.25). The differences between the 
models are very small for the highest T e g and lowest log g 
values. Otherwise they depend strongly on the wavelength 
range, and no general trend is visible. This is illustrated 
by Table 2, which lists the (T e g, log g) combinations for 
three metallicities in order of increasing flux differences 
from top to bottom. Three different wavelength ranges 
have been regarded: blue, UV and red. It can be seen 
that in the latter two, the convectivc efficiency effect is 
inversed compared to the one in the blue part. Thus, one 
can only guess that the calculated emitted flux depends in 
a complex way on the combination of T e fj, logg, and the 
convection model. 

We use our grids of computed fluxes to derive colors 
and color indices in the uvby photometric system. The role 
of convection on this photometric system has already been 
studied by SK (see Sect. 2), also using the ATLAS9 code 
in the MLT and CM cases. Here, we extend this study 
to the CGM case, and investigate how the variations of 
color indices due to temperature and gravity variations 
arc affected by the convection formulation. 

We find that: 

— there are no measurable differences between colors or 
indices computed with CM and CGM models and both 
are very close to those computed with MLT(a=0.50) 
models. 



2 Information on the solar activity level has been ob- 
tained from the National Solar Observatory Digital Library 
(http: / /www.nso. noao.edu/diglib /ftp. html) . 
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Fig. 8. Top: Observations of the solar irradiance from Thuillier et al. (1998), Lockwood et al. (1992) and Neckel & Labs (1984) 
(solid lines) and solar irradiance calculated with the CGM model (dashed line) . Bottom: Solar irradiance calculated with three 
different convection models. 



— differences become much more important if colors 
or indices computed with MLT(a=1.25) models are 
compared to those computed with CM, CGM, or 
MLT(a=0.50) models. 

Thus, we can extend the results of SK (see Sect. 2) 
to the CGM model, and conclude that color indices 
computed with CGM models are generally in better 
agreement with observations than those computed with 
MLT(a=1.25) models. 

The (b — y) index is the most sensitive one with re- 
spect to temperature changes, and this sensitivity is also 
strongly influenced by the convection model considered. 
We have investigated the variation of (b — y) indices com- 
puted using models differing only by the convection treat- 
ment. Fig. 10 shows that the sensitivity of (b — y) to 
convection change is T c g and gravity dependent, and the 
temperature changes associated with the ones of (b — y) 
are written along the curves. The results are very similar, 
and in the same order of magnitude, for mctallicities ten 
times or one tenth of the solar one. The same conclusion 
is reached when the CM model is replaced by the CGM 
or by MLT(a=0.50) formulation. 



From this result we can establish that the "er- 
ror" (or "change") on temperature variation estimations 
can be important, i.e. as large as 200 K, when using 
MLT(er=1.25) instead of CM convection treatment (or 
CGM or MLT(a=0.50)). It can reach up to 400 K, if the 
overshooting option of ATLAS9 is not removed. 



6. Conclusions 

One of the main conclusions to be raised from this study is 
that as long as one considers inefficient convection, what- 
ever is the choice of the formulation, either MLT with 
low a, or FST, the interpretation of spectroscopic or pho- 
tometric observations is equivalent: observed BLPs and 
Stromgren color indices of dwarf and subgiant stars be- 
tween A5 and G5 spectral types, and in a large range 
of metallicity are best represented by the use of less effi- 
cient convection transport, i.e. MLT with a = 0.5, or with 
FST formulation. This confirms results already obtained 
by Fuhrmann et al. (1993), VM and van 't Veer-Mcnncrct 
et al. (1998) for the Sun, Procyon, and other cool metal- 
poor stars using MLT models. 



U. Heiter et al.: New grids of model atmospheres I 17 




300 400 500 600 700 

I [nm] 



Fig. 9. Ratios of MLT and CGM to CM fluxes for three dif- 
ferent combinations of T e tt and log<?. [M/H]=— 1 for all cases. 
Relative transmissivity profiles of the Stromgren uvby filters 
are indicated (multiplied by a standard 1P21 photomultiplier 
response function). 



Table 2. Model parameters ordered according to increasing 
flux differences arising when using any two different convection 
models. The quantities listed in the columns labeled "blue" , 
"UV" and "red" are the maxima of |-ffA,MLTi.25/#A,CM — 1| 
for the wavelength ranges 360-520 nm, 250-360 nm and 520- 
700 nm, respectively (cf. Fig. 9). 
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Fig. 10. Differences between (b — y) indices computed using 
models differing only by the convection treatment. The thin 
line is for log g = 2.5, the thick line for logg = 4.5. For clarity, 
we indicate the temperature differences corresponding to the 
(b — y) differences along the curves only for a few models. 



Gardiner et al. (1999) reported a few opposite cases 
(see Sect. 2), but for parts of their sample of stars funda- 
mentally known log<? values were not available. An anal- 
ysis of a larger sample of stars in binary systems with 
revised fundamental parameters for both \ogg and T e g 
(Smalley et al. 2002) did not confirm the discrepancies pre- 
viously found. Furthermore, for the case of F stars Smalley 
et al. (2002) noticed a larger systematic difference between 
fundamental effective temperatures and those obtained 
from H^j lines for MLT (a= 1.25) than for less efficient con- 



vection models, although this discrepancy remains within 
the overall uncertainties. 

Nevertheless, we have to emphasize that in mod- 
els with deep convection zones (e.g. for Sun, Procyon) 
MLT(cv = 0.5) and FST treatments have comparable ef- 
fects on calculated fluxes, but not on atmosphere struc- 
ture. They produce different temperature gradients in 
the deep layers, as can be seen in Figs. 2 and 3, but 
those cannot be distinguished by the computed BLPs. In 
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other words, the BLPs allow to discriminate among dif- 
ferent values for the MLT parameter a, but not among 
MLT(a = 0.5), CM, and CGM models. In any case, the 
sensitivity of BLPs to convection parameters depends sig- 
nificantly on the other physical parameters. This holds 
especially for the sensitivity to gravity change, which can 
be more important than usually expected. 

In case of weakly efficient convection, fluxes and colors 
depend only weakly on the selected convection treatment. 
On the other hand, when the convection is highly efficient, 
then fluxes and colors become strongly dependent on the 
convection modelling, as the differences among the models 
show up more clearly within the photosphere. Thus, signif- 
icant uncertainties on stellar global parameters arise from 
the convective treatment in model atmospheres. Ignoring 
these uncertainties can lead to systematic differences af- 
fecting subsequent interpretations. 

The calculations of color and limb darkening par- 
tial derivatives are significantly improved when using the 
present model atmosphere grids which are finer spaced in 
T e ff and logg and have a higher resolution in the tem- 
perature distribution with depth (Barban et al. 2002). 
Smoothness of these derivatives is of crucial importance in 
the mode discrimination problem for non-radially pulsat- 
ing stars, which is basically due to the dependence of the 
color amplitude ratios on these derivatives. Details of the 
required precision of these partial derivatives in order to 
be useful for mode identification will be given in Garrido 
et al. (2002). There we will show that the next space as- 
teroseismology missions - COROT, MONS/R0mer and 
Eddington - will supply light curves with high enough 
precision to permit a direct comparison up to the second 
order to partial derivatives with respect to temperature 
and gravity as calculated with the present model atmo- 
spheres. 

The improved resolution of the new model grids also 
avoids unphysical oscillations in evolutionary track calcu- 
lations when using ATLAS9 model atmospheres as bound- 
ary conditions (see Sect. 4.1). Moreover, the possibility to 
choose among different convection models allows a self- 
consistent match between model atmospheres and model 
envelopes (Montalban et al. 2001). However, we must 
stress here that the different relations T and F conv /F to t 
vs. depth represent stars which are different in their radii 
and luminosities. The broad effects of the convective treat- 
ment can only be assessed by studying a complete stellar 
model, i.e. a model with an atmosphere and an internal 
structure which are consistently built with the same con- 
vection formulation. We will address this topic in follow 
up work (Kupka et al. 2002). 
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Appendix: Description of SMGT 

The program can be run in either of two modes depending 
on the temperature structure used for initialization: 

- In the static mode, an existing model file or a gray 
atmosphere is used. 

— In the dynamic mode, an existing model file or a 
weighted average of all existing models within one grid 
step of each parameter is used. The weights are calcu- 
lated in the following way: For each atmospheric pa- 
rameter p, the quantity 

r_jL;£L\ (27) 

\ Pma,x Pmin / 

is computed, where i denotes the initialization model, 
m the model to be computed, and p ma x/min the max- 
imum and minimum of the parameter as given in the 
grid definition. The results for all parameters (at most 
four) are multiplied to yield the weight. 

The state of convergence of a particular model is mea- 
sured by calculating the RMS and the maximum of the 
flux errors (AF/), the flux derivative errors (AF/) and the 
temperature correction (AT;) in all layers I. A model is 
considered as fully converged if these values satisfy certain 
criteria, which are given in Table 3 (labeled "primary"). 
Models for which these criteria cannot be achieved after 
a reasonable number of iterations are also stored, if they 
satisfy the criteria labeled "secondary" in Table 3, but the 
corresponding files are marked with "~" . 

In order to achieve the convergence criteria without 
wasting time when no further improvements can be ex- 
pected from further iterations, the required number of it- 
erations is calculated and checked dynamically, after an 
initial sequence of 12 iterations. From a sequence of n 
iterations the ATLAS9 output is processed every n/4 it- 
eration (but at least every 15 th and at most every 3 rd 
iteration) and the speed of convergence is characterized in 
the following way. The ratios between the rms errors of 
two subsequent iterations (r F ,r F ,) are calculated. If they 
are found smaller than a threshold value (0.95), damping 
exponents are computed iteratively for the flux errors: 

J F =[7 F - 1 n i - 1 -aln(l-r i F )]/n\ (28) 
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where i goes from 2 to the number of processed outputs, 
jp is set to zero or the value determined from the previous 
iteration sequence, a = 1 for i — 2 and a = [n l — ri L ~ l )/ri l 
for i > 2, and an analogous formula is used for the flux 
derivative errors. The total number of iterations is then 
given by 

Inherit) 



n tot = int 



or ntot 



int 



ln(l — cxp(— 7i?n)) 

ln(r crit ) 
ln(r>) 



with 



^crit 



RMS prim (AF,) 



(29) 



(30) 



RMS(AF ; ) 

where the largest of the values resulting from AFi or AF( 
is taken. Apart from the first iteration sequence, the actual 
number of further iterations is a fraction depending on 
the ratio between the current predicted number of total 
iterations and the previous prediction. The computations 
are terminated if the errors increase (after a few further 
trials) or if the total number of iterations would be too 
large (we use a limit of 480). 
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